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I. INTRODUCTION 



The Delta baryons, with spin and isospin 3/2, have long been a subject of study 



y, s q, y, y, 0. 



n. 



M in QCD Sum Rules (QCDSR), which are a well-established method 

nn 

that links the hadron degrees of freedom with the underlying QCD parameters [9|, [Hj. An 
important issue in the analysis of the decuplet baryons as well as of other excited hadrons 
is their finite widths. Conventionally, the method is to apply the zero-width approximation, 
where the Breit-Wigner contributions to the imaginary parts of the correlation functions 
are represented by 5-functions 0, [h]]. However, one cannot automatically assume that 
the contributions of the finite-width effects are negligible. Therefore, it is necessary to 
go beyond analytical arguments and investigate numerically the finite-width effects on the 
hadron spectrum in QCDSR. 

The finite-width contributions within QCDSR have been considered before in the con- 
text of the p- meson 11] and the p-uo mixing 12, [l3]. Our aim in this work is to make a 
reanalysis of the Delta sum rules by taking into the account the finite-width effects. For this 
purpose, we first consider the 'traditional' sum rules and in order to implement the finite- 
width effects on the A(1232) state, we use the method developed in Ref. [yj], which was 
then applied to the first pion excitation state, 7r(1300). Our second aim in this paper is to 
include the first positive- and negative-parity Delta excitation states, A(1600) and A(1700) 
respectively, within QCDSR analysis as sub-continuum resonances and to extract informa- 
tion about the parameters that characterize these resonances. For this purpose we resort 
to the parity-projected sum rules in order to separate the positive- and the negative-parity 
contributions [3, Q, Q- 

For our numerical procedure, we use the Monte-Carlo based analysis introduced in 
Ref. jl7j ]. This method provides a more systematic treatment of uncertainties in QCDSR. 
On the other hand, we expect that a Monte Carlo based analysis can better resolve the 
discrepancies in the hadron spectrum due to finite-width contributions. We have organized 
our paper as follows: in Sec. [Ill we present the formulation of the traditional Delta sum 
rules including the finite-width contributions. In Sec. IHH we consider the parity-projected 
sum rules and compare the results with those from the traditional ones. Finally, we arrive 
at our conclusions in Sec. HVl 
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II. TRADITIONAL SUM RULES 



Construction of the Traditional Sum Rules 



We begin our analysis by considering the correlation function 



n M „ = % I d*x e ip - x ( 



(1) 



where rj is the Delta interpolating field: 



Vi*( x ) = e abc [u a r (x)C'j ll u b (x)}u c (x). 



(2) 



Here a, b, c are the color indices, T denotes transposition and C = ry 2 7°. The QCDSR are 
derived by calculating the correlation function in Eq. ([T]) in the deep Euclidian region using 
Operator Product Expansion (OPE) and matching this to a phenomenological ansatz. The 
correlation function in Eq. ([!]) has the following form: 



IV (p) = Tli{p 2 )g^ + n 2 (p 2 )^ + . . . 



(3) 



where the ellipsis represents the Lorentz-Dirac structures other than and g^j). In prin- 
ciple, one can obtain the sum rules at different Lorentz-Dirac structures. Here, we are 
interested in the g^ v and the g^ft structures, which are contributed by only the J = 3/2 
particles (see e.g. Refs. [l, 5] for a more explicit separation of the spin components). 

The analyticity of the correlation function allows us to write the phenomenological side 
of the sum rules in terms of a dispersion relation of the form 

i r imn 

7T 



Ren^(g 2 



■ ds. 



(4) 



The ground-state hadron contribution is singled out by utilizing the zero-width approxima- 
tion, where the hadronic contributions from the Breit-Wigner form to the imaginary part of 
the correlation function is proportional to the 5-function: 

-1 



[- 



lim Im . 

r^o Ls — m 2 + imT 



ttS(s — m ), 



Imn^(s) =n'£ t 6(8-m 2 x )(0\r tll \X)(X\T h ,\0). 



(5) 
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One then expresses the correlation function in Eq. ([I]) as a sharp resonance plus a continuum 
after Borel transformation: 

1 r°° 

ni (M 2 ) = \lme- m2 / M2 + - lmU ie - s / M2 ds, 

1 \1 W2 (6) 
n 2 (M 2 ) = ~\le- m2 ' M2 + - / Imn 2 e- s / M2 d S) 

7T J w 2 

at the structures and g^ft, respectively. Here, M is the Borel mass, w is the continuum 
threshold and the matrix element of the current ff between the vacuum and the Delta state 
is defined as 

(0|^|A(p,s)> = W(p,s), (7) 

where Aa is the overlap amplitude and v^ip, s) is the Rarita-Schwinger spin vector of the 
Delta. We also make use of the Rarita-Schwinger spin sum, which is 

J2 «"(p, sW(p, s) = - (sT - \fi> - ^f/^ - 2 ~^f) U + m). (8) 

s 

For the sake of completeness, we first present the sum rules in momentum space: 

3 3 18p 2 J (p 2 - s) 

(W) : £ ln( V) -S'W^.^^ < < = *[^) *■ 
Upon Borel transformation, the chiral-odd sum rule (OSR) at the structure g^ u is given as 

-aExL^M 4 - ^E m 2 aL 2 / 27 M 2 - —abL 16 / 27 ^ A 2 m / I(m, T) e~ s ^ M2 ds, (10) 
> 3 18 Jo 

and the chiral-even sum rule (ESR) at the structure g^j) for the Delta is given as 



(9) 
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-E 2 L^ 27 M 6 - —b E L^ 27 M 2 + -na 2 L 2 *' 27 - -^m 2 a 2 L 17 ' 27 
5 72 3 9M 2 

f w 2 (11) 

= A 2 / /(m, T) e- s/A/2 ds, 
io 

where A = (27t) 2 Aa . To include the finite-width effects, we replace the 5-functions 

in Eq. (jSj) by a Breit-Wigner form fis| 

7 ( m > H = -7 2^2, 2r2 ; (12) 

7r (s — m z ) z + m^l 1 

where T is the width of the resonance state. In the sum rules above we have defined the quark 
condensate a = —(2Tr) 2 (qq), the gluon condensate b = (g 2 G 2 ), and the quark-gluon mixed 



condensate (qg c a ■ Gq) = m^qq) with the QCD coupling-constant squared g 2 = 4ira s . 
The four-quark condensate is parameterized as ({qq) 2 ) = n(qq) 2 - The corrections that 
come from the anomalous dimensions of various operators are included with the factors 
L = log(M 2 /AQ CD )/log(/x 2 /Ag CD ), where /i = 500 MeV is the renormalization scale and 
Aqcd is the QCD scale parameter. The perturbative corrections are taken into account with 



the factors 
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dimension-five (m a): 1 



dimension-zero (/): 1 + ^ 
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(14) 



in ESR, where 'Je — 0.58 is the Euler constant. These corrections give large contributions 
to ESR (leading-order correction is ~ 70% for a s /ir ~ 0.12 at the scale of 1 GeV 2 ) whereas 
their contributions to OSR is small and therefore can be safely neglected. The continuum 
contributions are represented by the factors E n {x) = 1 — e~ x 'Yl n x n /nl, with x = w 2 /M 2 . 

B. Numerical Analysis of the Traditional Sum Rules 



In order to obtain an analytic approximation for the integral over the Breit-Wigner shape, 
we first express the imaginarypart of the Breit-Wigner shape as a Riemann sum of unit-area 
rectangular pulses P m (s,T) [141 ] : 



P m (s, T) = [B(s - m 2 + mT) - Q(s - m 2 - mT)]/2mT, 



mT 



lim — 



(s — m 2 ) 2 + m 2 T 2 n^oo n ^—^ y j — f 



1 Pm[S, 



n 



J-f 



ir, o</<i, 



(15) 



where an approximation for the integral over Breit-Wigner shape can be obtained by using 
(for w 2 > m? + m r) 



P m (s,T)e- s ' M2 ds 



2 /M 2 



G(m,T,M 2 



M 2 / 
G(m, T, M 2 ) ee — sinh 
mT \ 



"mT 
M 2 



(16) 
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TABLE I: The obtained values of the parameters representing A (1232) from a consideration of 
1000 parameter sets and using OSR in Eq. ([10]) and ESR in Eq. (fTT]) . The values of the search 
parameters are given both in the zero-width limit (T — > 0) and with finite width (T = 0.12 GeV). 
The fixed parameters are denoted by an asterisk in each case and the numbers inside brackets are 



experimental values as given by PDG 



20] 



Region (GeV) 


T (GeV) 


m (GeV) 


w (GeV) 


A 2 (GeV 6 ) 




[0.118 ± 0.002] 


[1.232] 






ESR 


* 


1.21 ±0.26 


1.65 ±0.25 


1.69 ± 1.04 


1.15 < M < 1.35 


0.12* 


1.24 ±0.21 


1.60 ±0.23 


1.61 ±0.83 




0.28 ±0.09 


1.232 * 


1.57 ±0.08 


1.34 ±0.29 




0.27 ±0.09 


1.23 ±0.12 * 


1.55 ±0.13 


1.36 ±0.37 


OSR 





1.43 ±0.12 


1.65 ±0.22 


2.39 ± 0.94 


0.95 < M < 1.10 


0.12 * 


1.48 ±0.11 


1.72 ±0.21 


2.80 ± 0.99 




0.16 ±0.06 


1.48 * 


1.70 ±0.04 


2.69 ±0.39 




0.20 ±0.08 


1.48 ±0.15 * 


1.76 ±0.20 


3.12 ± 1.18 



In order to proceed, we have made a 4-pulse approximation by taking n = 4 in Eq. (|15|) and 
choosing / = 0.7 in order to make the area of the four pulses equal to the total area under 
the Breit-Wigner curve (= tt) 1 . One can then write the integral over the Breit-Wigner shape 
as 



- f 7 ^~ 2r2 e~ s/M2 ds = e- 2 /^ w[ r M 2 ]; (17) 



where 



W[m, T, M 2 ] = 0.1592 G(m, 3.5119 T, M 2 ) + G(m, 1.4412 T, M 2 ^ 



(18) 



+ G(m, 0.8597 T, M 2 ) + G(m, 0.4606 T, M 2 ) 

Note that W — > 1 in the narrow-resonance limit (r — > 0). This four-pulse approximation 
serves to reduce the numerical difficulties in evaluating the integrals over the Breit-Wigner 
shape, which extends to low-s region. We will further discuss this issue in Section IIHI 

We determine the uncertainties in the extracted parameters via the Monte Carlo based 
analysls Educed in Re, In this a nalysls , ran d omly selecte d, Ga.sia.ly distributed 



1 In the n — > oo limit, Eq. (fTSJ) would be true for any choice of / between and 1 
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sets are generated from the uncertainties in the QCD input parameters. Here we use a = 
0.52 ±0.05 GeV 3 , b = 1.2 ±0.6 GeV 4 , m 2 = 0.72 ±0.08 GeV 2 , and A QCD = 0.15 ±0.04 GeV. 
The factorization violation in the four-quark operator is searched via the parameter k, where 
we take k — 2 ± 1 and 1 < k < 4; here (qq 2 ) > (qq) 2 is assumed via the cut-off at 1 (for 
a discussion on QCD parameters see e.g. Ref. [l7J). We use 1000 such configurations from 
which the uncertainty estimates in the extracted parameters are obtained. 

We begin our analysis by making a three-parameter fit to the left-hand sides (LHS) of 
the sum rules in Eq. ( TTUj) and Eq. (fTTj) . including A, m and w. The valid Borel regions are 
determined so that the highest-dimensional operator contributes no more than 10% to the 
OPE side in order to warrant the OPE convergence, while the continuum contribution is less 
than 50% of the phenomenological side, which provides pole domination. Note that, while 
the first criterion is rather straightforward, one does not initially have a complete control 
on the second, since the phenomenological parameters are determined from the fit and they 
are correlated. We use the following strategy: we first make the fits in a reasonably selected 
Borel region, which is then adjusted by trial and error according to the fit results, until the 
above criteria are satisfied. We also seek a region which best resolves the resonances and 
finally obtain an optimized Borel region. 

In Table HJ we present the obtained values of these parameters from a consideration of 
1000 parameter sets. The values of the search parameters are given both in the zero-width 
limit (r — > 0) and with finite width (r — > 0.12 GeV) [20]. In the zero-width approximation, 



the Delta mass values we obtain are in agreement with those in Ref. [6J. It is well-known 

that OSR performs better as compared to ESR even though OSR somewhat overestimates 

the Delta mass. We observe from the fitted parameter values in Tabled that the finite-width 

effects lead to a change of ~ 3% in the Delta mass whereas these effects change the continuum 

thresholds by 3-5%. Moreover, the uncertainties in the values of the fit parameters from ESR 

are reduced by ~ 20% with the finite-width effects. We have also made a three-parameter 

fit including T, A and w with the input m = 1.232 GeV for ESR, and m = 1.48 GeV for 

OSR. From ESR, we obtain an overestimated value for the width as T = 0.28 ± 0.09 GeV, 

which should be compared with the experimental result T = 0.118 ± 0.002 GeV. On the 

other hand, OSR successfully produces the experimental width value. Note that the overlap 

amplitude values we obtain from the Monte-Carlo fits are in agreement with the ones from 

Jl 

earlier determinations: A 2 = 2.3 ±0.6 GeV 6 from the traditional analysis of the sum rules [2|], 
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X 2 = 2.0 



model 



22 



0.7 GeV 6 from lattice QCD [2l| and A 2 = 1.8 ± 0.5 GeV 6 from instanton liquid 



The main conclusion of Table [U is that the mass determination is rather insensitive on 
the width variation. The masses in the zero-width approximation and with the finite-width 
effects agree within their uncertainties. At this point, it is also important to analyze the 
correlations between the phenomenological parameters, in particular between the mass and 
the width. The truncated OPE does not have enough information to allow a simultaneous 
determination of these two parameters. In order to extract one of them, we require the other 
as an input from experiment. However, given the approximate nature and limited accuracy 
of QCDSR (considering the truncated OPE, uncertainties in the condensate values and the 
continuum model), one needs to consider variations in all input parameters, including the 
mass and the width. For this purpose, we have allowed a variation of 10% in the mass values 
by generating Gaussianly distributed sets from Monte Carlo, and studied the correlations 
by taking these sets as input. Our results are given in Tabled we observe that a variation 
in the mass leads to larger uncertainties and deviations for the width in OSR, while this 
behavior is not evident for ESR. 



III. PARITY-PROJECTED SUM RULES 



A. Construction of the Sum Rules 



Our next task in this analysis is to include the first Delta excitation on the phenomenolog- 
ical side by assuming that both the lowest-lying and the first-excitation Delta states reside 
below the continuum threshold and all the other subsequent excitations can be embedded in 
the continuum. One of the obstacles with this consideration is the identification of the first 



excitation, since the interpolating fie 
also to the negative-parity state 15 



d in Eq. (j2J) couples not only to the positive-parity but 



16] . If one is interested in the lowest-lying positive- 
parity resonance, the excited states can be regarded as parts of the continuum. However, 
when the first excitation is taken as a sub-continuum resonance and included as a pole on the 
phenomenological side, the question arises whether this pole actually represents the positive- 
parity A(1600) resonance or the negative-parity A(1700) resonance. The complexity of the 
excited-state spectrum may not allow us to separate the resonances of interest from the 
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continuum by a proper choice of the threshold. Therefore, in the following, we apply parity 
projection to the correlation function in Eq. (Tj[|) in order to separate the positive- and the 
negative-parity contributions [3, Q, Q so that we can include the first excitation in our 
QCDSR analysis. For this purpose, we use the old-fashioned correlation function, which is 



II*, = i J d A x e ip ^(xo)(o|T[^(x)^(0)]|o). (19) 



In the zero-width resonance approximation, one can write the imaginary part of the 
correlation function in the rest frame p = as (for the g^ and the g^p 1 structures) 



Im rP(po) = £ [(A+) 2 ^-5(Po - m+) + (AJ 2 ^-5(p - m r , 

n 

= l0 A{p ) + B{p ), 
where A(p ) and B(p ) are defined as 

A (Po) = \ £[(A+) 2 <KPo - m+) + (\-) 2 5(p - m-)}, 

n 

B(Po) = ~ £[(A+) 2 <HPo ~ O - (KWpo ~ m-)], 



(20) 



(21) 



n 



and A± = (27r) 2 A±. The parity- projected Delta sum rules are then given as 
A(M, w±) ± B(M, it>±) = 

fw± rw± 

-P 2 /M 2 



2m ± (A ± ) 2 / I(m±, r ± ) dpQ + 2m > ± ( y ±) 2 / J(m ^ r ^ 

Jo Jo 



)e- p 3/^^ Pl 



(22) 



o- 



with 



A(M, w ± ) = 1 P± L 4 / 27 - A 6 P± L 24 / 27 + | « a 2 L 28 / 27 , 
B(M, w ± ) = A -a Pt L 16 / 27 - 2 - ml a P± Z 2 / 27 , 



(23) 



where we define 

P±(M,w±)= I dp p n e-Po/ M \ (24) 



w± 



Note that the functions A and .B can be readily identified from Eq. ([9]). In Eq. (1221) . m + 
(m_) and A + (A_) denote the mass and the overlap amplitude of the positive- (negative- 
parity Delta resonance, respectively, w + is the continuum threshold for the positive- 
(negative-)parity sum rule and the primed notation refers to those parameters for the first 
positive- (negative-)parity excitation. We have also made use of the relation 

S(p - m±) = 2m± 5 (pi - m|) (25) 



on the positive po-axis and replaced the 5-function with the Breit-Wigner form in Eq. (TT2]) 
(taking s —>■ p^). We have included the first Delta excitation with the second term on the 
right-hand side of the sum rule in Eq. (1221) . The functions in Eq. (1231) can be written more 
explicitly as 



A(M,w ± ) = - 



144 
2 

B(M,w±) = - a 



b (l - e- w ^ M2 ) M 2 L 24 / 27 + | K a 2 L 28 / 27 , 



L 4/27 



(26) 



M 2 L 16 / 27 -^m 2 a/(^ ± )L 2 / 27 , 



with I(w) = f™ e~ x2 l M2 dx. 



B. Numerical Analysis of the Parity-Projected Sum Rules 

We follow a similar procedure as above and first try a three-parameter fit to the LHS 
of Eq. (122]) for the positive-parity state including m + , A + and w + , with the zero-width 
approximation and assuming that the first Delta excitation lies in the continuum. The 
valid Borel regions are estimated similarly as in the traditional sum rules: we take 0.95 < 
M < 1.30 GeV. A three-parameter fit returns a continuum threshold value smaller than the 
resonance mass, which is an unphysical solution. This case, which signals that the OPE does 
not have enough information to resolve all the three parameters simultaneously, has been 
realized before in various other works on QCDSR with Monte-Carlo fits 12J, [2|. In 
order to proceed then, we switch to a two-parameter fit by fixing the continuum threshold at 
a value suggested by the the Particle Data Group (PDG) [20(] and the traditional sum rules 
above, as was also done in the previous works. To put it more explicitly, before fixing it 
we have an idea about the value of the continuum threshold from the traditional sum-rules 
analysis and experimental Delta spectrum. We fix the continuum threshold at a reasonable 
value accordingly. 

At this point we would like to make a technical remark: the physical spectrum of the 
Delta resonance with a p-wave starts at po = (m^ + m n ) while the tail of the infinite Breit- 
Wigner shape we utilize extends unphysically to low-po region. In order to eliminate the 
overlap of the spectrum with the region that falls below the n-N threshold, one can use 



a parametrization of the spectral function in terms of a pcrdependent width (see Ref. 11] 
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for the corresponding case of the p-meson). Such a constraint aggravates the numerical 
difficulties that arise from a multi-parameter fit. Therefore, we choose to use the simplest 
parametrization of the Breit-Wigner form with a constant width 20j, which serves best for 
our purposes. For the parity-projected sum rules, the part of the spectrum below the tt-N 
threshold amounts to 15% of the total spectrum (for physical mass and the width). In the 
case of the traditional sum rules, where the physical spectrum starts at s = (itln + m 7T ) 2 , we 
have smaller overlap thanks to the four-pulse approximation in Eq. ffToT) : actually for OSR, 
which performs better as compared to ESR, all the rectangular pulses are above the n-N 
threshold, while for ESR, only 9% of the Breit-Wigner distribution lies below the threshold. 

In Table |TTJ the obtained values of the fit parameters from a consideration of 1000 param- 
eter sets are presented, for two different values of the continuum threshold; w + = 1.5 GeV 
and w + = 1.6 GeV. In the zero- width approximation for w + = 1.5 GeV, we obtain 
m + = 1.13 ± 0.09 GeV, whereas a two-parameter fit including m + and A + with the finite- 
width effects gives m + = 1.23 ± 0.09 GeV. For the latter, we evaluate the integrals over 
the Breit-Wigner shape numerically and obtain a Delta mass value in agreement with the 
experimental one. If we instead fix the continuum threshold at w + = 1.6 GeV, the mass 
values produced by the fit are changed by 10%. Comparison of the mass values obtained 
with the zero-width approximation and with the finite-width effects suggests that the finite- 
width contributions are rather important for the parity-projected Delta sum rules. These 
effects lead to a change of ~ 20% and ~ 10% in and m+, respectively. A search for 
the width and the overlap amplitude in a fit with m + = 1.232 GeV successfully reproduces 
the experimental result; that is T + = 0.118 ± 0.002 GeV 20J]. The overlap amplitudes are 
compatible with those in Table [I] from the traditional sum rules. In order to analyze the 
impact of the mass variation on the width, we allow a 10% variation of the input ground 
and excited state masses. We observe that the uncertainties in the widths are increased 
by a mass variation. To further illustrate the correlation of the mass and the width of the 
A(1232) resonance, in Fig.(CQ) we present the results obtained with a broad range of mass 
and width using OSR and the parity-projected sum rule. As was shown in Section [Til OSR 
performs better as compared to ESR, however overestimates the mass value. It is seen from 
the figure that the masses start to deviate from their initial values towards higher width 
region for both sum rules. It is clear that the mass increases with the increasing width 
and the lowest mass is given by the calculation of the spectral function with a 5-function 
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FIG. 1: Correlation of the mass and the width of the A(1232) resonance using OSR (filled squares) 
and the parity-projected sum rules (PPSR) (empty squares). The diamond marks the experimental 
mass and the width of A(1232). Each datum point has been obtained from a consideration of 100 
configurations. The error bars associated with each point stem from the uncertainties in the 
condensate values. 



representation. This observation is in agreement with that of Ref. ll|] in the case of the 
p-meson. 

In order to include the first positive-parity excitation in our analysis, we apply a three- 
parameter fit by taking m + = 1.232 GeV, T + = 0.12 GeV and m' + = 1.6 GeV. In this 
case, the upper limit on the valid Borel window is increased according to the selection 
criterion above, since the continuum contributes less with the increasing threshold: we take 
0.95 < M < 1.50 GeV. In Table [Til we present the width values of the first positive- 
parity excitation obtained from the fit, for different values of the resonance mass and the 
continuum threshold. We also allow a reasonable variation in the input mass and continuum 
threshold parameters, by generating them via Monte Carlo. The resulting width values are 
in agreement with the experimental result, which is V, = 0.35 ± 0.10 GeV as given by 



PDG 20[. A comparison of the overlap amplitude values produced by the fit indicates that 
the lowest-lying and the first excited states couple to the interpolating current with strengths 
of the same order. 
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TABLE II: The obtained values of the parameters representing the positive-parity A(1232) and 
A(1600) resonances from a consideration of 1000 parameter sets using the parity-projected sum 
rule in Eq. (|22|) . For the excited A (1600) resonance, the mass and the width of the lowest-lying 
A(1232) resonance are fixed at m + = 1.232 GeV and r + = 0.12 GeV, respectively. The Borel 
window is taken as 0.95 < M < 1.30 GeV when we consider only the A(1232) resonance, while 
the window is broadened to 0.95 < M < 1.50 GeV when the A(1600) resonance is included. The 
fixed parameters are denoted by an asterisk in each case and the numbers inside brackets are 
experimental values as given by PDG 2C|. 



A(1232) 


T+ (GeV) 


m + (GeV) 


w+ (GeV) 


A^GeV 6 ) 




[0.118 ±0.002] 


[1.232] 








* 


1.13 ±0.09 


1.5 * 


1.26 ±0.19 




* 


1.19 ±0.09 


1.6 * 


1.57 ±0.20 




0.12 * 


1.23 ±0.09 


1.5 * 


1.48 ±0.19 




0.12 * 


1.31 ±0.09 


1.6 * 


1.83 ±0.23 




0.14 ±0.08 


1.232 * 


1.5 * 


1.46 ±0.28 




0.10 ±0.06 


1.232 * 


1.6 * 


1.98 ±0.47 




0.15 ±0.10 


1.23 ±0.12 * 


1.5 * 


1.55 ±0.49 


A(1600) 


r'+ (GeV) 


m' + (GeV) 


w' + (GeV) 


(A' + ) 2 (GeV 6 ) 




[0.35 ±0.10] 


[1.625 ± 0.075] 








0.22 ±0.07 


1.6 * 


1.8 * 


3.47 ± 0.45 




0.20 ±0.09 


1.6 ±0.1 * 


1.8 * 


3.14 ±0.87 




0.20 ± 0.08 


1.6 * 


1.85 ±0.05 * 


3.42 ± 0.73 



In this work, we have also analyzed the negative-parity Delta states; however, the scarcity 
of the experimental information for the negative-parity resonances prevents us from con- 
straining the parameters, which represent these resonances. Particle Data Group (PDG) 
lists the lowest-lying negative-parity (J = 3/2) state as A(1700) with Breit-Wigner width 
of 0.3 ± 0.1 GeV. Assuming that the first Delta excitation lies in the continuum, we have 
first applied a two-parameter fit to the LHS of Eq. ( 1221) including m_ and A_ and have 
found that the experimental mass value of the A(1700) resonance can be obtained only for 
w_ > 2.4 GeV (see Table HIT]) . Fixing the continuum threshold value at u>_ = 2.4 GeV, 
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TABLE III: Same as Table U but for the negative-parity A(1700) and A(1940) resonances. For 
the excited A(1940) resonance, the mass and the width of the lowest-lying A(1700) resonance 
are fixed at m_ = 1.7 GeV and T_ = 0.3 GeV, respectively. The Borel window is taken as 
1.30 < M < 1.55 GeV. 



A(1700) 


T_ (GeV) 


m_ (GeV) 


w- (GeV 6 ) 






[0.30 ±0.10] 


[1.710 ± 0.040] 








* 


1.45 ±0.67 


2.4 * 


0.68 ±0.41 




0.30 * 


1.57 ±0.32 




0.95 ± 0.29 




0.42 ± 0.36 


1.70 * 




0.95 ±0.18 


A(1940) 


r'_ (GeV) 


m'_ (GeV) 


w'_ (GeV 6 ) 


(A'_) 2 




[?] 


[« 1.940] 








0.65 ±0.29 


1.94 * 


2.4 * 


1.41 ±0.88 



we obtain m_ = 1.45 ± 0.67 GeV with the zero-width approximation. If we adopt the av- 
erage value for the width of the A(1700) resonance as T_ = 0.3 GeV [20] and apply a fit 
with finite-width effects, we obtain m_ = 1.57 ± 0.32 GeV, while a search for the width 
with m_ = 1.7 GeV produces T_ = 0.42 ± 0.36 GeV. The relatively high value, at which 
we have to fix the continuum threshold in order to obtain the experimental mass value of 
the A(1700) resonance suggests that the excited-state contributions may not be embedded 
in the continuum and therefore the QCDSR analysis cannot be limited to the lowest-lying 
resonance. The first negative-parity excitation (with J = 3/2) is listed as A(1940) by PDG 
and the available data for the Breit-Wigner width of this resonance have large uncertain- 
ties (0.10 ~ 0.78 GeV). To continue, we include in our fit the A(1940) resonance with 
m'_ = 1.94 GeV. We apply a three-parameter fit by fixing the continuum threshold value 
at w~ = 2.4 GeV and by taking m_ = 1.7 GeV, T_ = 0.3 GeV. As a result of this fit we 
obtain V'_ = 0.65 ± 0.29 GeV, which favors a large value for the width of A(1940). Despite 
the naive analysis as a result of the large uncertainties in the phenomenological input, our 
results for the negative-parity Delta resonance parameters are consistent with those given 
by PDG |2o|. 
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IV. CONCLUSION AND DISCUSSIONS 



In summary, we have analyzed the finite-width effects on the QCD sum rules of the 
Delta baryons, through explicit utilization of the Breit-Wigner shape. We have observed 
that although the finite-width effects give minor contributions to the traditional sum rules 
(about 3% change in the A(1232) mass value), the effects on the parity-projected sum 
rules are rather significant (about 10% change in the A(1232) mass value). The deviations 
from the narrow-resonance approximation in QCDSR analysis leads to results in better 
agreement with the experiment for the parity-projected case. Our findings for the width of 
the A(1232) resonance from a Monte Carlo analysis of the parity-projected sum rule are in 
agreement with the experimental results. The first positive-parity Delta excitation, A(1600), 
has also been considered as a sub-continuum resonance using parity-projected QCDSR, and 
its Breit-Wigner width has been calculated, where the experimental width value has been 
produced. For both resonances, in order to analyze the mass-width correlation, we allowed 
a variation in the input mass and the width values. We have found that while the masses 
are insensitive on a width variation for small width values, the uncertainties in the widths 
are increased with respect to a mass variation. This in turn means that the boundaries are 
more restrictive for the mass and less for the width. Our numerical analysis shows that the 
masses increase with the increasing widths and the lowest mass is given by the calculation 
of the spectral function with a 5-function representation. For the negative-parity case, large 
experimental uncertainties prevent us from constraining the input parameters, therefore only 
a naive analysis can be made. A reasonable agreement with the experimental values has 
been obtained though, which favors a large value for the Breit-Wigner width of the A(1940) 
resonance. It is desirable to extend this analysis to all decuplet baryons and their excited 
states. A work along this line is still in progress. 
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